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Most of our effort has been devoted to preparing a paper on 

the dissociation cross section of OH arising from absorption into 
2 - 

the repulsive 1 T, state. A copy of the paper, THE LOW-LYING 
2 - 

E STATES OF OH, by Ewine F. van Dishoeck, Stephen R. Langhoff 

and A. Dalgarno, is attached. It has been submitted for publication 

to the Journal of Chemical Physics. 

We have also examined again photodissociation by predissociation 
2 + 

of the A E state. A pi eliminary version of a paper describing 
our results is attached as part of this report. 
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Abstract 


2 - 

Potential energy curves for the two lowest I states of OH are computed 

at the configuration-interaction level using four different basis sets. Elec- 

tronic transition dipole moments connecting the excited 1 E and 2 (D) E 

2 

states with each other and with the ground X n state are presented as 

functions of intemuclear distance. The theoretical absorption oscillator 

2 - 2 

strengths for the D S (v'=OJ X H(v"=0) transition are in good agreement 
with the empirical value derived from astronomical measva-ements . The photo- 
dissociation cross sections for absorption from the v"=0, 1 and 2 levels of 

2 “ 

the ground state into the continuum of the 1 E state are calculated, and 
the interstellar and cometary photodissociation rates are derived. 
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1. Introduction 

The hydroxyl radical OK is an inportant component in a diverse array of 

1 

physical environments. We have reported recently theoretical calculations 

2 + 

of the properties of the excited E electronic states of OH in which the 

2 + 2 + 2 + 

oscillator strengths and radiative lifetimes of the A E » B E and C S states 

were derived. We present now calculations at a similar level of accuracy 

2 - 2 - 

for the two lowest S states. The lowest 1 E state is known from predisso- 
• 2 + . 

ciation of the A E state. The state is repulsive and may provide an import- 
ant path for the dissociation of OH in cometary atmospheres and in the 

interstellar gas. There are no measurements of the absorption cross sections. 

2 ** 2 “ o 

The second 2 E state is the spectroscopic D E state in the 1200 A region 

2 

studied in absorption by Douglas . Observations of the absorption lines 

have yielded estimates of the oscillator strength and the abundance of OH 

. 3 

in diffuse interstellar clouds . 

Calculations of limited accuracy of the potential energy curves of the 

l^E and D^E states have been carried out by Lefebvre-Brion et al.^ and 

5 2 - 2 4 

Easson and Pryce . The D E state was classified as a Rydberg state ' . No 

2 - 2 

calculation has been reported of the 1 E -X H transition moment. The 
2-2 

D E -X n transition moment has been computed at one Anternuclear distance 

g 

by Ray and Kelly in a ten-configuration self-consistent field calculation. 

2. Theoretical Methods 

We employed the configuration-interaction (Cl) method for the. determina- 

2 - 

tion of the electronic wavefunctions of the E states using four different 
atomic orbital basis sets. Two of the basis sets, denoted I and II, consisting 
of Slater functions, and a third, denoted III, consisting of Gaussian 
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2 - 

functions, are given explicitly in Ref. 1. Since the D S state contains 
substantial Rydberg character, an adequate representation requires diffuse 
basis set functions, in particular 3s and 3po functions centered on the 
oxygen nucleus. Basis sets II and III satisfy this requirement, but basis 
set I, which contains a diffuse s function only on the hydrogen nucleus, does 
not. Accordingly, we constructed a fourth basis set, denoted lA, that 
augments basis set I by adding diffuse 3s (exponents ^ = 1.0 and 0.5), 4s 
(5 = 0.5 and 0.2), 3pa (? - 1.0 and 0.5) and 3pn (? = 0.75) functions cen- 
tered on the oxygen nucleus. The Is (C = 0.4) and 2s (5 = 1.557) functions 
centered on the hydrogen nucleus were removed, to avoid linear dependencies. 

The resulting 22a llir 86 basis should represent accurately both the compact 
and diffuse components of the wave function. 

In the calculations with basis sets I and lA, the internally-consistent 

7 2 

self-consistent field (ICSCF) X II molecular orbitals were used for all the 
states. The configuration state functions (CSFs) that were explicitly in- 
cluded in the wave function were selected by Rayleigh-Schrodinger perturbation 
theory from the list of CSFs generated as single or double excitations from a 
zeroth-order wave function, . Any CSF having a coefficient exceeding 0.08 

in the wave functions for any of the states at any internuclear distance was 

2 

included in ijj^. For the E states, the representation of contained 

2 + 2 - 

numerous terms, in part because the E and E states were considered together. 
In the calculations, those CSFs contributing more than 20 microhartrees (yh) to 
the energy were retained, and the energy contribution of those discarded 
wai estimated by an extrapolation procedure^. Although these methods provide 
an accurate baleinced description of the excited states, an uncertainty of 
about 0.05 eV exists in the energies resulting from the extrapolation pro- 
cedures, and from difficulties in obtaining a uniform zeroth-order 
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representation cf the states as the internuclear distance varies. 

2 - 

In all Cl calculations with basxs set II, we used the 1 S molecular 

2 + 

orbitals as in the computations on the S states. These orbitals have the 

2 

advantage over the X H state orbitals that an additional orbital (4o) is 

2 

stabilized by the Hartree-Fock method. For the ground X II state and the 
2 - 

first excited 1 S state, a first-order Cl (FCXJl) calculation was performed which 
incorporated between 300 and 700 configurations^. For the FOCI calculation the mol- 
ecular orbitals were partitioned into core orbitals la, internal orbitals 20-4o 

2 - 

and Itt, and external oibitals 5a-24a, 2h-10tt and 16-26. For the D Z state, the 
internal set consisted of the 2a-10a and liT-4Tr orbitals and the external set 
contained the orbitals lla-24a, Stt-IOtt and 16-26. A set of reference con- 
figurations was constructed which included all configurations with 
Cl coefficients greater than 0.1 at any internuclear distance. All 
single excitations from th'e reference configurations and 
double excitations within the internal set were retained but double excita- 
tions out of the 2a cjrbital were excluded. The selection process resulted in a 
Cl matrix with 2400 configurations. The energies of this more elaborate calcula- 
tion were shifted to the FOCI energies so that both yielded the-saiae excitation 

2 - 

energies for the 1 Z state. 

2 - 

The Cl calculations with basis set III also used the 1 Z orbitals, and 
employed the Wuppertal-Bonn MRD-CI package of programs for the construction 
of the. Cl matrices. These programs select the CSFs on the basis of their 
energy lowerings in a secular equatpLon with the reference configurations but the 
procedure is otherwise analogous to that described previously for bases I and lA. 
The selection threshold was set at 50 )ih. An improved energy estimate was ob- 
tained by correcting the extrapolated Cl energies for the effect of unlinked 
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8 9 2 + 

clusters ' . This led to negligible improvement for the Z states and was 

not done in the earlier calculations^. Additional details of our Cl methods are 

given in Ref. 1. 

The bound -bound and bound-free absorption cross sections for transitions 
between electronic states depend on the electronic transition moment D (R) . If 
the z-axis of the electron coordinates lies along the internuclear axis, D(R) 
is the expectation value, in atomic units, of the operators Ez. or E (x ,+iy . )//2" 

j ^ j ^ ^ 

for transitions in which |AA| = 0 or 1, respectively. In terms of D(R) , the 

2 

electronic absorption oscillator strength at R for transitions from the X II 
2 - 

State to the E states is given by 

f^^(R) = I gAE(R) |d(R)|^ (1) 

where g is a weighting factor which equals unity for E transitions 

and Ae(R) is the energy difference at R in hartrees. The band 

absorption oscillator strength for the transition from a bound vibrational 

2 - 

level (v") of the ground state to a bound vibrational level (v* ) of a E 
state is defined by 


v"v' 


gAE, 


v'v" ’ V' v' 


(2) 


where AE^,^„ is the transition energy measured in hartrees and is 

the transition matrix element 


^V'V" (X^||D(R) Xyll) 


(3) 
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obtained from the normalized vibrational wave functions Xy(R) the initial 
and final states. The vibrational wave functions were determined by 

nvimerical integration of the radial Schrodinger equation for nuclear motion 
ignoring centrifugal distortion. 

The absorption cross section from a bound vibrational level v” of the 

2 - 

ground state to the vibrational continuum of a S upper state can be written 

o^„ = 1. 23x10”^ g|(Xj^(R) |d(R) |x^„(R)) |^AE cm^ (4) 

where k is the wave number of the continuum state. In terms of k, the kinetic 

2 

energy of relative motion in the final state is Ej^ = k /2y where y is the 
reduced mass, all quantities being in atomic units. The continuum wave 
function normalized to the asymptotic form 

1 

Xk<^) ^ ^|k^^ sin(kR+Ti) , (5) 

where n is an elastic phase shift. The spontaneous transition probability for 

the emission of a photon of wavelength X % from a bound upper state vibra- 

2 - 

tional level v' to a lower Z continuum state with wavenumber k" is given by 

20 2 

|(Xj^„(R)|d(R)|x^, (R))| (6) 

measured in photons s 
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3, Results and Discussion 
3.1 Potential Energy Curves 

2 2 - 2 - 

The relative potential energies of the X H, (1) Z and D Z states of 

OH calculated with basis sets I and lA, are presented in Table I and those 

corresponding to basis set i are illustrated in Fig. 1. The character- 

2 

istics of the ground state X H potential resulting from calculations I-III 

1 2 - 

have been discussed in detail previously . The 1 E potential is repulsive 

everywhere. The principal spatial occupancy occurring in the description of 

2 2 2 — 

this state is (lo) (2o) (3a) (4o)(lir) (lit) at all internuclear distances, and it 

separates correctly at infinite distances into ground state oxygen and 

hydrogen atoms. In the separation, the 4a orbital changes from an antibonding 

OH orbital at a distance of 1.85 a^ to a hydrogen Is orbital at infinity. 

2 

The calculated vertical excitation energy from the minimum in the X It state 

2 “ 

(R„ = 1.8342 a ) to the 1 Z state is 8.00, 7.97, 7.20 and 7.81 eV for 
o 

basis sets I, lA, II and III wave functions, respectively. The good agreement 

between the basis set I and lA results indicates that Rydberg functions are 

2 - 

unimportant to the description of the 1 E state. The basis set I and lA 

2 

energies tend to be high because they use the X II orbitals, whereas the 

2 - 

basis set II and III results may be low both because they used the 1 E orbitals 

and because smaller Cl wave functions were employed. From the larger Cl cal- 

* 2 - 

culations, it appears certain that the vertical excitation energy to the 1 E state 

2 - 

is 7.95 +0.1 eV. The precise location of the 1 E state is crucial to the 

determination of the OH photodissociation rate in cometary atmospheres. 

2 - , . 

The D E state is described principally by the spatial cccupancy 
2 2 2 

da) (2a) (3a) (5a)(lTr) (Itt) at 1.85 a^, where the 5a orbital contains significant 

3s and 3pa Rydberg character. At larger distances, other spatial occupancies 
2 2 — 

such as (la) (2a) (3a) (4a) (5a) (lir) (lit) become important, in part because the 


9 


molecular orbitals are themselves nOt optimum for the description of the 

2 “ 2 “ 3 

D state. The D S state dissociates to the excited 0( S) state of 

oxygen and a ground sta-te hydrogen atom. In contrast to the 
2 - 

repulsive 1 2 state, it is strongly bound. Using the experimental oxy- 

10 2-2 
gen atom energy levels and the excitation energy T^ of the D 2 state and dis- 

2 IX 

sociation energy D of the X n state , we obtain a binding energy of 3.95 eV 
6 

2 “ 

for the D E state. The calculated spectroscopic parameters T^ and for 
2 - 

the D 2 state are compared with experiment in Table II. The basis set I 

2 

results are in remarkable, though in part fortuitous, agreement with experiment , 

whereas the basis set lA results for both T^ and are slightly too small. 

The basis set II wavef unctions overestimate R and underestimate T , whereas 

e e 

basis set III wave functions underestimate both, particularly R . The good 

agreement obtained for T^ with basis set I arises from a cancellation of 

errors. On the one hand, basis set I does not contain sufficiently diffuse 

2 - 

3s and 3pa functions to describe accurately the D 2 scate and on the other 

2 

hand it omits a larger amount of the correlation energy of the X n state 
2 - 

than of the D 2 state. Our methods yield about the same percentage, 70% for 

2 2 - 

basis set I, of the valence correlation energy of the X n and D 2 states, 

2 

but becav^se the total correlation energy is greater for the compact X H 

2 - 2 - 
state than for the diffuse D 2 state, a bias is created in favor of the D 2 

state. When the basis set deficiencies are removed by using basis lA, the 

predicted values of T^ and R^ are too small. This problem is more severe 

for the basis set III calculations, which account for the smaller fraction, 

60%, of the valence correlation energy and for the basis set II calculations, 

where only 20% of the valence correlation energy was recovered. Similar 

6 

arguments apply to the limited MCSCF calculations of Ray and Kelly , 
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3.2 Transition Moments 

2-2 2-2 

The calculated electronic transition moments of the 1 E -X n and D E -X II 

transitions obtained with basis sets 1, lA, II and III are given in Tablet,’. Ill 

and IV and are illustrated in Figs. 2 and 3, respectively. The basis set I 
2 - 2 

results for the 1 E -X II transition are not shown in Fig, 2. They are almost 

identical to the basis set lA results, demcjnstrating the insensitivity of the 

transition moment function to the addition of diffuse functions. The close 

agreement between all of the theoretical transition moment functions, especially 

for R values near the ground state minimum, indicates that the magnitude of the 

theoretical photodissociation cross sections should be reliable. 

2 - 2 

The theoretical D E -X H transition moment functions in Fig. 3 show 

larger discrepancies. There are substantial differences between the basis set 

I •’\nd lA results, particularly for R < 2.25 a^ where the diffuse character 
2 - 

of the D E state is most prevalent and reduces the magnitude of the transition 
moment. The basis set lA transition moment function is in reasonable agree- 
ment with the basis set II and III results for small internuclear distances, 

but rises more steeply at larger distances. Beyond 2.5 a the theoretical 

2 - 2 

transition moment functions stay approximately constant. The D E -x H moment 

should approach 0. 40 t 0,04 at infinite separation, the value appropriate to the 

3 3 X2 6 

0( P)-0( S) transition . The MCSCP value of Ray and Kelly computed at the 

equilibrium separation of the ground state is also shown in Pig. 3. It agrees 

well with the basis set II and III results. 

2 - 2 - 

The electronic dipole moments for the D S -1 K transition obtained 
with basis sets lA and II are given in Table V and illustrated in Fig. 4. 

The magnitude of the transition moment decreases rapidly with increasing R. 

The two theoretical transition moments are in good agreement for R £ 2.5 a^ 
where the transition moment has large amplitude. 
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3.3 Oscillatoj, strengths and photodissociation cross sections 

2 - 

a) Absorption into the D 2 state 

The electronic absorption oscillator strengths at the equilibrium 

separation of the ground state, and the 0-0 band absorption oscillator 

2 - 2 

strength for the D 2 -X n transition are tabulated in Table VI. The 
theoretical energy differences were used in computing f^^^. To determine the 
band absorption oscillator strength f^^, we used the calculated potential 
energy curves in the equations for the vibrational wave functions, together 
with the corresponding transition moment function. The potential energy 
curves were shifted to be consistent with the experimental excitation energy. 

The basis set lA transition moments were used in conjunction with the | 

ti 

basis set I potential curves. An empirical value of f^^ has been obtained | 

3 

from astronomical observations of OH in diffuse interstellar clouds . The de- p 

rived value of 0.0108 + 0.003 agrees closely with the theoretical values obtained 

with wave functions lA, II and III. The basis set I result of 0.016 is 

probably a reliable upper bound, because the lack of 3s and 3po oxygen basis 

2-2 

functions tends to produce an overestimate of the D 2 -X H transition 

13 

moment. These results rule out the earlier estimate of 0.1. 

Because the dipole moment changes rapidly with R, the close agreement 
obtained with wavefunctions II and III is partly fortuitous. If the basis 
set I potential energy curves, which are positioned closest to the experi- 
mental curves, are used together with the basis set II and III transition 
moments, changes of up to 30% in f^^ are observed. The rapid variation of 

D(R) is the major source of the inaccuracy in the prediction of Ray and j 

6 ' “ 
Kelly , which was carried out in advance of the astronomical determination. f 

2-2 

Table VII lists the absorption oscillator strengths f^^, of the D 2 -X II I 

I 

I' 
1; 

I 

Ij 

ii 


X 
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transition for upper state vibrational levels up to v'=3, calculated using the 

basis set I potential curves and basis set lA transition moment function. 

2 - 2 - 
The D £ state can radiate to the lower repulsive 1 £ state ^ which leads 

to dissociation of the molecule. The continuous spectra emitted in the tran- 
2 — 2 ~ 

sition from the D Z state to the 1 Z state are shown in Figure 5 for v'=0 

and 1, obtained using the basis set I potential curves and the basis set lA 

transition moment function. The integrated transition probability, is 

7-1 7-1 

3.4x10 s for v'=0 and 3.2x10 s for v'=l, values which agree to within 

10% with those obtained with the basis set II potential curves and transition 
2 - 

moments. The D £ state can also radiate into the ground electronic state. 

2 

The summed transition probability into the discrete levels of the ground X II 

8 —1 8 "X 

state we calculate to be 2.3x10 s for v'=0 and 2.4x10 s for v'=l using 

the basis set I curves and the basis set lA moments. With the basis set II 

8-1 8-1 

results we obtain 1.5x10 s for v'=0 and 1.6x10 s for v'=l. Thus, about 

2 - 

one fifth of the absorptions into the D S v*=0 and 1 levels lead to disso- 
ciation of the molecule. With the oscillator strengths and wavelengths of 
Table VII; the dissociation rate corresponding to the interstellar radiation 
field used by Roberge at al; ^^is 1.0x10 s ^ in the optically thin limit. 

2 - 

b) Absorption into the 1 Z state 

2 - 

Absorption into the 1 E state leads directly to dissociation. The 

• -4 15 

effective oscillator strength has been estimated as 8x10 by Solomon and 

•2 16 

as 10 by Tarafdar and Das. We obtain at R = 1.85 a a value of 

o 

—2 2 — 

3.4x10 for f^^ so that absorption into the 1 Z state is a major destruction 

mode for OH in diffuse interstellar clouds and in comets. 


13 


The photodissociation cross sections obtained with wave functions lA, II 

2 

and III for absorption from the v"=0 level of the X II state are illustrated in 
fig. 6 as a function of incident photon energy. The cross section shapes ob- 
tained with the three wave functions are similar but, consistent with the dif- 
ferent vertical excitation energies, there are small shifts in the effective 
thresholds. Table VIII contains the numerical values obtained with basis 
set lA. 

The cross sections listed in table VIII lead to a photodissociation rate 
""6 “"1 

of 2.2 X 10 s at the mean Earth-Sun distance (1 A.U.) in. conditions of average 
solar intensities, the rate changing by less than a factor of two between 
solar maximum and minimum. The solar spectrum is varying rapidly in the 

2 - 

vicinity of the effective threshold for absorption of photons into the 1 2 
state. The uncertainty of 0.1 eV in the calculated excitation energy pro- 
duces an uncertainty of 25% in the solar induced photodissociation rates. 

The interstellar radiation field is nearly uniform in the spectral 

2 - 

region where absorption into the 1 Z state occurs, and the interstellar 
photodissociation rate is insensitive to the detailed location of the 

14 

threshold. Corresponding to the radiation field used by Roberge et al., 
we obtain in the optically thin limit a rate of 1.3 x 10 s 

17 

In experimental situations, OH is often vibrationally excited. 

Pigs. 7 and 8 illustrate the photodissociation cross sections for absorption 
from the excited v"=l and v"=?2 levels of the ground state using the basis 
set lA potential curves and transition moments. The structure reflects 
that of the initial vibrational wave function. 
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Table I 



The X^n, 1^ 

- 2 - 
S 2ind ■> 2 

potential 

energies 

in eV calculated 



using basis set 

I and lA 



R(a^) 


x^n 

i=: 

s” 


2 ~ 
D r 


I 

lA 

I 

lA 

I 

lA 

1.30 

4.07 

- 

- 

- 

- 

- 

1.50 

1.18 

1.20 

10.03 

10.00 

12.05 

11.76 

1.75 

0.06 

0.06 

8.31 

8.29 

10.57 

10.28 

1.8342 

0.00 

0.00 

8.00 

7.97 

10.37 

10. 10 

2.00 

0.14 

0.15 

7,55 

7.50 

10.22 

9.98 

2.25 

0.70 

0.74 

6.94 

6.91 

10.31 

10.12 

2.50 

1.37 

1.43 

6.35 

6.33 

10.63 

10.48 

2.75 

2.06 

2.12 

5.86 

5.82 

11.05 

10.97 

3.00 

2.66 

2.68 

5.46 

’ 5.47 

11.48 

11.46 

3.25 

3.16 


5.16 


11.96 


3.50 

3.54 


4.95 


12.40 


3.75 

3.82 


4.80 


12.77 


4.00 

4.03 


4.68 


13.10 


4.50 

4.26 


4.55 

■ 



5,00 

4.38 


4.48 




6.00 

4.47 


4.43 
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Table II 


2 - 

Spectroscopic parameters for the D S state of OH 


Parameter 

Basis 

Set 


and g 
Kelly^ 

Experiment 

I lA 

II 

III 

R (a_) 
e o 

2.05 2.02 

2.06 

1.95 


2.043 

Tg (eV) 

10.22 9.98 

9.16 

9.93 

<9.4 

10.20^ 


^The value of 10.20 eV was obtained by correcting the reported Tq value by 

2 - 2 

the difference in zero-point energies between the D E aund X II states. 
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Table III 


2 “ 2 a 

Theoretical 1 E -X II electronic transition moments of OH in a.u. 


R(ao) 



Basis Set 


I 

lA 

II 

Ill 

‘1.00 

- 


0.480 

- 

1.20 

- 

- 

- 

0.481 

1.30 

- 

- 

0.482 


1.50 

0.456 

0.463 

0.472 

0.468 

1.75 

0.422 

0.430 

0.439 

0.433 

1.8342 

- 

0.412 

- 

- 

1.85 

- 

- 

0.416 

- 

2.00 

0.353 

0.363 

0.368 

0.359 

2.25 

0.272 

0.278 

0.276 

0.254 

2.50 

0.209 

0.215 

0.205 

0.174 

2.75 

0.165 

0.170 

- 

- 

3.00 

0.131 

0.136 

0.126 

0.098 

3.25 

0.105 

- 

- 

- 

3.50 

0.084 

- 

0.080 

0.061 

3.75 

0.066 

- 

- 


4.00 

0.050 

- 

0.048 

0.036 

4.50 

0.029 

- 

- 

- 

5.00 • 

0.012 

- 

- 

- 

6.00 

0.003 



0.002 


^The operator is (x+iy)//2^ 
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Table IV 


2—2 a 

Theoretical D E -X H electronic transition moments for OH in a.u. 




Basis 

Set 


R(a^) 

I 

lA 

II 

III 

1.30 

- 

- 

. 0.127 


1.50 

0.218 

0.153 

0.148 

0.143 

1.60 

- 

- 

0.162 

- 

1.75 

0.276 

0.216 

- 

0.175 

1.8342 

- 

0.246 

- 

- 

1.85 

- 

- 

0.211 

0.197 

2.00 

0.351 

0.310 

0.247 

0.238 

2.25 

0.402 

0.394 

0.301 

0.308 

2.50 

0.412 

0.435 

0.329 

0.347 

2.75 

0.402 

0.446 

- 


3.00 

0.411 

0.456 

0.331 

0.354 

3.25 

0.384 


- 


3.50 

0.352 


0.332 

0.339 

3.75 

0.324 


- 

- 

4.00 

0.291 


0.340 

0.291 


Ray and Kelly 


6 


0.190^ 


^The operator is (x+iy)/»^ 

^Their result in the velocity formulation is Q.179 
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Table V 


Theoretical D^E 

2 - 

-IE electronic 

transition 

moments 

R(a^) 


Basis Set 


lA 


II 

1.30 

4 


2.373 

1.50 

2.292 


2.310 

1.60 

- 


2.239 

1.75 

2.018 


- 

1.8342 

1.831 


- 

1.85 

- 


1.886 

2.00 

1.399 


1.536 

2.25 

0.815 


0.887 

2.50 

0.501 


0.490 

2.75 

0.195 


- 

3.00 

-0.299 


0.141 

3.50 



-0.134 

4.00 



-0.218 

8.00 



-0.482 
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Table VI 

2 ~ 2 

The absorption oscillator strengths of the D Z -X H transition 




Basis 

Set 


g 

Ray and 

Observation^ 


1 

lA 

II 

III 

Kelly 

a 

el 

0.023 

0.015 

0.010 

0.010 

0.008 

- 

b 

00 

0.016 

0.012 

0.011 

0.011 

0.004° 

(0.0108 + 0.003) 


^See Eq. 1 for definition. The theoretical excitation energies were used. 

See Eq. 2 for definition. The experimental T^ value was used. In the 
vibrational analysis, the theoretical potentials were used. 

Obtained assuming a constant transition moment, and a Franck-Condon factor 
of 0.45, based on a Morse representation of the potentials. 
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Table VII 

2 “ 2 

Band eibsorption oscillator strengths and wavelengths for the D E (v')-X n(v"«0) 



transitions^ 

0 

V’ 

^Ov' 

A (A) 

0 

1.2(-2) 

1222.1 

1 

4.4 (-3) 

1188.7 

2 

5.8(-4) 

1158.4 

3 

9.6{“6) 

1130.9 


^Calculated using the basis set I potential curves and the basis set lA 
transition moment function. 
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Table VIII 


2 - 

Calculated photodissociation cross sectic^ns for absorption into the 1 2 state from 
2 

the X H(v''=0) level using the basis set lA potential curves. and transition moment function 


E(eV) 

a (10 cm 

6.50 

0.01 

7.00 

0.28 

7.20 

0.70 

7.40 

1.42 

7.60 

2.28 

7.80 

2.90 

7.90 

3.01 

8.00 

2.97 

8.20 

2.54 

8.40 

1.84 

8.60 

1.14 

8.80 

0.59 

9.00 

0.25 

9.50 

0.01 


Figure Captions 


2 2 - 2 - 

Fig. 1 - Theoretical basis set I potentials for the X 11/ 1 S and D E 
states of OH. 

2 “ 2 

Fig. 2 - Theoretical 1 E -x n electronic transition moments. The A's, d's and 
o's represent the basis set lA, II and III results/ respectively. 

Fig. 3 - Theoretical D E -X II electronic transition moments. The V'S/ A's, 

C's and o's represent the basis set I, lA, II and III results/ 

6 

respectively. The MCSCF result of Ray and Kelly at R = 1,8342 a^ in 
the length formulation is denoted by the x. 

2 - 2 — 

Fig. 4 - Theoretical D E -1 E electronic transition moments. The A's and o's 
represent the basis set lA and II results/ respectively. 

Fig. 5 - The continuous emission spectra originating from the v'=0 and v'=l 

2 - 2 - 
levels of the D E state in the transition to the 1 E state. 

Fig. 6 “ The photodissociation c*. >ss section for absorption from the v"=0 

2 2 - 

level of the X II state into the 1 E state calculated with wave func- 
tions lA/ II and III. 

Fig. 7 - The photodissociation cross sections for absorption from the v"=l lev- 

2 2 - 

el of the X II state into the 1 E state calculated with wave function 

lA. 

Fig. 8 - As for Pig. 7 for the v"=2 level. 
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The rates of photodissociation of the OH and OD molecules from 

2 2 + 

absorption of solar radiation in the X II-A E electronic transition 

are calculated to lie between 2.0 and 3,7x10 s for OH and to 
-7 -1 

be about 2.4x10 s for OD at 1 AU from the Sun. The corresponding 
lifetimes which are upper bounds to the actual lifetimes are general- 
ly consistent with the observational cometary data. 
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1. Introduction 

Jackson (1980) has pointed out that the OH radical in cometary 

atmospheres is destroyed by a photodissociation process initiated 

by the absorption of solar radiation into high vibrational levels of 

2 + 

the excited electronic A Z state. The vibrational levels of the 
A Z state with quantum numbers v' ^2 undergo predissociation with 
an efficiency close to unity for all rotational levels (Brzozowski, 
Erman and Lyyra 1978) . Because the dissociation is stimulated by 
line absorption in a wavelength region where the solar spectrum is 
rapidly varying the corresponding lifetime of OH depends upon the 
radial velocity of the comet (Jackson 1930) . The deuterated molecule 
OD is destroyed by a similar process at, as we shall argue, a much 
slower rate. We take the opportunity to present improved estimates 
of the photodissociation rate of OH. 

2. Photodissociation of OH 

Detailed calculations of the populations of the rotation- 

2 

vibration levels of the ground electronic X II state of OH (Mies 1974, 

Biraud et al. 1974, Elitzur 1981, Despois et al. 1982, Schleicher 

and A' Hearn 1982) show that at a distance of 1 AU from the sun, 

about 97% of the OH molecules reside in the J"-3/2, N"=l, v''=0 level 
2 

of the X HI 2/2 calculations are consistent with observa- 

tional data and in particular with high resolution ultraviolet 
spectra obtained from the lUE satellite (Jackson et al. 1979, 

Feldman et al. 1980) . 

For each band, six allowed absorption lines occur labelled 
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ORIGINAL PAG2 Sfi 
OF POOR QUALITY 


(1) , ^P 22 ^( 1 ) ^ Q|l^(l)/ ^Q 2 i( 1 )/ excitation 

rate for the absorption of solar radiation into vibrational level 
2 + 

v' of the A Z electronxc state is given for unit density of OH by 
the sum over the individual lines 


I Bq , W,N') P^,(|r s"^ 

J',N' ^ 


where p^, is the mean radiation density of the solar flux in units 

•2 

cm of the transition from the 


of ergs cm at the wavenumber v * in — 


level N"=l, v*' = 0) of the X^TI state to the level (J' , N' , v' ) 

2 + 3 

of the A S state and Bq^, l-J' , N' ) is the corresponding Einstein 


2 —1 —1 

B coefficient for absorption in units of cm ergs s . 


The B coef- 
3 


ficient may be written in terms of the coefficient A^,q(J', 1) 

for spontaneous emission according to 




In terms of the absorption oscillator strength fQ^i ( ) , 


n 7 2 




To a close approximation,’ the line absorption oscillator strength 
may be expressed in terms of the band oscillator strength hy 
the relationship 


S(J"N"-J'N') ^ 
(2J"+1) ^Ov' 


fov' (J”N”-J'N') 
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where S ( J"N”“J'N’ ) are Honl‘-London factors norinalized according to 
^ S(J"N"-J'N') _ , 

j„ (2J"+1) ~ ^ * 

In the expression used by Jackson (1980) for the Einstein B coef- 
ficient, the Honl-London factors appear to be normalized incorrectly. 

Many estimates have been made of the absorption oscillator 
2 + 2 

strengths of the A 2 -X II transition in OH. Boring, Wadt and 

Cartwright (1982) have reviewed the experimental and theoretical 

data and they recommend a radiative lifetime c,f. 688 ns for the 
2 + 

A Z state which corresponds to an oscillator strength fQQ of 
1,02x10 . Recently McDermid and Laudenslager (1982) have measured 

a lifetime of 721 + 9 ns corresponding to f^^^ = (0.976 + 0.012) x 10 
Their value may be affected by rotational excitation. Wang and 
Killinger (1979) and Chidsey and Crosley (1980) have shown that 
fjjQ decreases slowly with increasing rotational quantum number. For 
the rotationless limit Wang and Huang (1980) recommend the value 
(1.09+0.04) X 10 ^ and Smith and Crosley (1981) the value 
(1.10 + 0.03) X 10"^. 

The values of fQ^i for higher values of v' are uncertain. The 
limited experimental data (cf. Crosley and Lengel 1975) suggest that 
foi/foo about 0.3 but no reliable information exists on fQ 2 * 
Because of extreme cancellation in the evaluation of the vibrational 
matrix elements the theoretical predictions are sensitive to details of 
the transition dipole moment as a function of intemuclear distance. Table 1 lists 
the calculated band oscillator strengths corresponding to the dipole 
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moment of Chu et al. (1974) and to the dipole moments of Langhoff 

et al. (1982) obtained with their case I and case III wave functions. 

Comparison with the experimental data on the 0-0 transition suggests 

a slight preference for the values derived from the dipole moment 

of Chu et al. (1974) and we adopt them. 

-5 

The value of 5,4x10 for £ q 2 is equivalent to a spontaneous 

5 -1 

vibrational emission transition probability of 1.06x10 s 

The approximate estimate made earlier by Jackson (1980) is equiv- 

4 -1 

alent to the much smaller rate of 1.05x10 s . The discrepancy is 
attributable mainly to the variation with internuclear distance of 
the transition moment which is not taken into account in the pro- 
cedure followed by Jackson (1980). 

The Honl-London factors S(J"N”-J'N') depend upon the coupling 
2 

constant a of the X state (cf. Dxeke and Crosswhite 1962, 

Schleicher and A' Hearn 1982). For v” = 0, a = -7.547 (Dieke and 
Crosswhite 1962, Moore and Richards 1971). We present in table 2 
the Honl-London factors and the wavelengths of the six absorption 
lines from the J”=3/2, N"=l, v”=- level to the optically allowed 
J' , N' , v*=2 levels. 

The solar radiation densities at the transition wavelengths 
have been calculated by Jackson (1980) for a range of cometary 
radial velocities between + 60 km s”^. Using them we may cal- 
culate the rates of photodissociation of OH by absorption of solar 
photons into the v’=2 level. The values are compared with the 
results of Jackson (1980) in table 3. At a distance of 1 AU from 

the Sun, they range from 3.16x10 ^ s ^ for a radial velocity towards 

—X ”6 "*1 

the sun of 58 km s to 1.20x10 s for a radial velocity away 


“1 

from the sun of 60 km s . They are smaller than the rates cal- 
culated by Jackson (1980) by factors of about four. 

The contribution to photodissociation from absorption into the 
v'=3 level was considered briefly by Jackson (1980). The correspond 
ing absorption line wavelengths are listed in table 2. The ratio of 
the absorption oscillator strengths f 03/^02 about 0.15 and the 

ratio of the solar intensities near 245 nm and 261 nm is about 0.6 
(cf. Mount and Rottman 1981) so that the inclusion of the v'=3 
channel increases the photodissociation rates by 9%. Because the 
increase is small, we ignore any dependence on the radial velocity 
of the comet and simply add to the rates for v’=2 listed in Table 3 
a factor 0.09 of the zero velocity v'=2 value. Absorptions into 
levels v’>3 are negligible in comparison. 

3. Photodissociation of OD 

A similar analysis may be carried through for OD with the 

2 + 

interesting difference that the v'=2 level of the A E stats of OD 
does not undergo significant predissociation and it is absorptions 
into the v'=3 and v’-4 levels that determine the resulting photo- 
dissociation rate. We have calculated the absorption oscillator 

2 2 + 

strengths fQ^, of the X II-A Z transition of OD using empirical 
potential energy curves and the dipole transition moment function 
of Chu et al. (1974) . The values are given in table 4. 

Bergeman, Erman, Haratym and Larsson (1981) have carried out 
measurements of the lifetimes of rotational levels of the v'=0-3 
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vibrational states of OD. Their analysis of the data suggests a 

predissociation rate for low-lying rotational levels of the v'=3 

state of about 6.9x10^ s’*'^ out of a total rate of about 2x10 s 

The corresponding predissociation efficiency n is 0.34. The 

absorption lines occur at wavelengths near 2575 S (Clyne, Coxon 

and Woon Fat 1973) and the corresponding photodissociation rate at 

-7-1 

zero velocity is about 2.0x10 s . No estimate is available for the 
predissociation efficiency n for v'=4. We adopt a value of unity. 
The absorption lines lie near 2450 8. With an oscillator strength 
of 1.2x10 , the absorption into the v'=4 level leads to a photo- 
dissociation rate of approximately 3.5xl0"® The total photodissocia- 

tion rate of OD at 1 AU from the Sun through absorption into the 
2 + —7 -1 

excited A E state is 2.4x10 s , leading to the possibility of a 

substantial enhancement in the relative abundance of OD/OH compared 

to [D]/tH]. However there are other dissociation mechanisms, one 

2 - 

such being that due to absorptions into the repulsive 1 Z state 
(van Dishoeck, Langhoff and Dalgarno 1983), which will limit the enhancement. 

Because of the additional paths , the predicted lif e~ 
times from absorption into the A S state are upper bounds to the 
actual lifetimes. The analyses of cometary data on OH in comets 
(Keller and Lillie 1974, Blamont and Festou 1974, B'eldman and 
Brune 1976, Feldman et al. 1980, Festou 1981, Despois et al. 1981, 
Weaver et al. 1981) yielded lifetimes of the same order but 
generally below and hence consistent with the calculated values. 

A study of the additional destruction processes is in progress. 
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Table 1 


Theoretical band oscillator strengths* of OH corresponding 
to different calculations of the transition dipole moment. 

Langhoff et al. 


O 

< 

Chu 

Case I 

Case III 

0 

l.lxlO"^ 

1.2xl0”^ 

1.15x10"^ 

1 

S-Ixio""^ 

3.1x10'"'^ 

2.9x10"'^ 

2 

5.4x10"^ 

5.2x10“^ 

4.4x10“^ 

3 

8.5x10"^ 

8.5x10“^ 

6.6xl0“® 

4 

1.4xl0"^ 

1.5x10“® 

1.4xl0“® 
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Table 2 

Honl-London factors and transition wavelengths for 
X^n 3/2 (J"=3/2, N”=l, v"=0) " N’, v'=2, 3) 


Line 

Wavelength 

nm 

v'=2 

S (J"N"-J'N' ) 
(2J"+1) 

Wavelength 

nm 

v'=3 

p^d) 

261.63 

0.1471 

245.00 


261.42 

0.0979 

244.83 

Q^d) 

261.42 

0.1404 

244.83 

®Q2i(1) 

261.00 

0.0616 

244.49 

R^d) 

261.00 

0,0413 

244,49 

CO 

to 

260.38 

0.0117 

243.97 
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Table 3 

Photodissociation rates by predissociation of OH at 1 AU as a 


function of heliocentric radial velocity v in units of 10 



V (km s 
-58 
-46 
-31.9 
-23 
0 

46.6 


(v'=2) 

3.16 

3.46 

1.95 

2.58 

2.58 

3.24 

1.80 


Rate (10“^ s"^) 
total 
3.39 
3.69 
2. 18 
2.81 
2.81 
3.47 
2.03 


Jackson (1980) 
13.09 
14.58 
8.70 
11.36 
4.85 
13. 51 


60 


7.41 


V' 

^Ov’ 

*The 

of 
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Table 4 

Theoretical band oscillator strengths fQ^i of OD* 

0 12 3 4 

l.lxlO"^ 3.2x10“"^ 5.6x10“^ 8.2xl0”^ 1.2xl0“® 

/ 

V 

dipole moment of Chu et al. (1974) was used in the calculation 
Ov' * 


